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In the Horndeski's most general scalar-tensor theories with second-order field equations, we derive 
the conditions for the avoidance of ghosts and Laplacian instabilities associated with scalar, tensor, 
and vector perturbations in the presence of two perfect fiuids on the flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) background. Our general results are useful for the construction of the- 
oretically consistent models of dark energy. We apply our formulas to extended Galileon models in 
■ which a tracker solution with an equation of state smaller than —1 is present. We clarify the al- 

(N ■ 

lowed parameter space in which the ghosts and Laplacian instabilities are absent and we numerically 
confirm that such models are indeed cosmologically viable. 
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I. INTRODUCTION 



The problem of dark energy is one of the most serious problems faced by cosmologists and particle physicists. The 
q'' cosmological constant is a simplest candidate, but it is plagued by the huge difference of energy scales between the 
Q^, theoretical and observed values Many alternative models to the cosmological constant have been proposed to 

approach the origin of dark energy f5|. 
bJOl, Dark energy models can be broadly classified into two classes: (i) modified matter models, and (ii) modified gravity 
models. In the class (i) the modified matter source is introduced to drive the late-time cosmic acceleration. The 
I ' representative model in this class is quintessence in which the potential energy V{(f)) of a slow- rolling scalar field 
^ I (/) is the source for dark energy. Another model is k-essence [4], in which a non-linear term with respect to the field 
00 ' kinetic energy X = — (9(/))^/2 leads to the cosmic acceleration. 

: The representative models that belong to the class (ii) are those based on f{R) theories IE, theories 01, 

00 Brans-Dicke theories [8], Dvali-Gabadazde-Porrati (DGP) braneworld [9], and Galileon gravity [lOl. (see Ref. 

' for reviews). The Lagrangian in f{R) theories is an arbitrary function of the Ricci scalar R, which corresponds to 
the particular class of Brans-Dicke theories (Brans-Dicke parameter wbd = The f{R,Q) theories consist of a 

function of R and the Gauss-Bonnet term Q, in which case there are two independent scalar degrees of freedom in a 
general background fli"]. In Galileon gravity, the Lagrangian is constructed such that the field equations are invariant 
under the Galilean symmetry df^cj) d^<j)-\-hfj^ in the limit of Minkowski spacetime [10] . One of those terms is given by 
X □(/), which appears in the DGP model as a result of the mixing of a brane-bending mode with a transverse graviton 
[m . In the Galileon model there exists a stable late-time de Sitter solution where 4> = constant [l^ ■ Moreover 
^ . , the model can be consistent with local gravity constraints through the Vainshtein mechanism ^1^ . 

■ Recently Deffayet et al. [l^ derived the action for the most general scalar-tensor theories having second-order field 
equations. This action is equivalent to that derived by Horndeski [23| in 1974 in the context of Lovelock gravity [2l| 
(see also ReL [22J). In fact the Horndeski's action is sufficiently general to accommodate most of the dark energy 
models proposed in literature. Moreover, the Gauss-Bonnet couplings F{(j))Q [l^, the generalized Galileon term 
G{(l),X)U(j) [2i - [27t . the derivative couphngs C^"" d^(l)d„(j) with the Einstein tensor G^^" [Hi also belong to the class of 
the Horndeski's action [2l|. 

For such general theories the full background and perturbation equations were recently derived in Ref. [29] in the 
presence of non-relativistic matter. In the context of infiation several authors computed the power spectra of scalar 
and tensor perturbations (2l| as well as the non- linear parameter of primordial non-Gaussianities (sol - fs^ . 

If we modify gravity from General Relativity we need to worry about the presence of ghosts as well as Laplacian 
instabilities. In the presence of two perfect fiuids (radiation and non-relativistic matter) we shall derive conditions 
for the avoidance of ghosts and Laplacian instabilities associated with scalar, tensor, and vector perturbations in 
the Horndeski's theories. This is of particular importance for the construction of theoretically consistent modified 
gravitational models of dark energy. We apply our results to extended Galileon models in which a tracker solution is 
present. This tracker corresponds to the generalization of that found in Ref. [16] for the covariant Galileon. Unlike 
the covariant Galileon the equation of state of dark energy wde can take the values close to —1 during the evolution 
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from the matter era to the accelerated epoch. This property should allow for the tracker to be compatible with 
observations. 

This paper is organized as follows. In Sec. the background equations of motion are derived on the flat FLRW 
spacetime in the presence of two perfect fluids. In Sec. IIIII we find conditions for the absence of ghosts and Laplacian 
instabilities by deriving the second-order action for scalar, tensor, and vector perturbations. In Sec. II VI we apply our 
general formulas to kinetically driven dark energy models, which cover the covariant Galileon as a specific case. Not 
only we identify a theoretically consistent parameter space but also we numerically integrate the field equations to 
confirm the analytic estimation. Sec. |V]is devoted to conclusions. 



II. THE BACKGROUND EQUATIONS OF MOTION 

In order to avoid the Ostrogradski instability [s^ it is desirable to keep the equations of motion at second order in 
derivatives. The most general 4-dimensional scalar-tensor theories having second-order field equations are described 
by the Lagrangian [l^ 



c = J2c., (1) 



where 

£2 = Ki(t>,X), (2) 

£3 = -G3(<^,X)n0, (3) 

£4 = G4(0, X)R + G4,x [(□<^)' - (V^V,0) (V^W)] , (4) 

£5 = G5(0, X) G^, (V^W) - i G^.x [{n<jyf - 3(00) (V^V,0) (V^V» -f 2(V'^V„0) (V"V^(/)) (V^V^0)] . (5) 

Here K and Gi {i = 3,4,5) are functions with respect to a scalar field (p and its kinetic energy X = —d^(l)df^(j)/2, R 
is the Ricci scalar, and G^u is the Einstein tensor. Gi^x and G^.^ {i = 3, 4, 5) correspond to the partial derivatives of 
Gi with respect to X and respectively, i.e. Gi^x = dGi/dX and Gi^^ = dGi/dcj). The above Lagrangian was first 
discovered by Horndeski in a different form |20|], which is equivalent to the Lagrangian ([T]) ^21i] . 

In addition to the scalar field ((> we take into account two barotropic perfect fluids, whose energy densities are pA 
and pB- Then the total action is given by 

S = J d^x^{C + CA + CB) , (6) 

where g is the determinant of the metric g^^^, , Ca and £ b are the Lagrangians of two perfect fluids respectively. 

Let us consider a fiat FLRW background with the metric ds^ = —N'^{t)dt^ + a^(t)dx^, where t is the cosmic time, 
N{t) is the lapse function, and a{t) is the scale factor. Varying the action ^ with respect to the N{t) and a{t) 
respectively, we obtain 

2XK^x -K + &X^HG3,x ~ 2XG3,0 - QH^Ga + 2AH^X{Gi,x + XG^.xx) - UHX^G^^^x ~ GH^G^^^ 

+2H^X4> (5G5,x + 2XG5,xx) - 6H^X (3G5,0 + 2XG5,^x) = -pA - pb , (7) 

K - 2X(G3,0 + 4>G3,x) + 2{3H^ + 2H)Gi - l2H^XGi^x - AHXGi^x - SHXGi^x - SHXXGa^xx 

+2(4> + 2i70)G4,0 + 4XG4^00 + 4X{4>- 2H<i>)G4.4,x - 2X{2H^ + 2HH4> + 3H'^4>)G5^x - AH^X'^4> G5,xx 

+mX{X - i/X)G5,0x + 2[2{HX + HX) + 2,H^X]G^,^ + iHX^Gs,^^ = -pA - Pb , (8) 

where a dot represents a derivative with respect to t, H = d/a is the Hubble parameter, and pa,Pb are the pressures 
of the two perfect fluids. Varying the action ^ with respect to (j){t), it follows that 

li(a3.)=P„ (<,) 



where 



J = ^K^x + QHXG3^x -2<j>G3,^ + 6H^iG4,x +^XG4,xx) -l2HXG4,4,x 

+2H''X{3G5,x + 2XG5^xx) - GH^iG^,^ + XG^,^x) , (10) 

= X0-2x(g3,00-I-0G3>x) +6(2if2 + i/)G4_^ + 6i/(i: + 2i/X)G4,0x 

-6if2xG5,00 + 2H^X^ G5,0x . (11) 
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The two perfect fluids obey the fohowing continuity equations 



PA 
Pb 



- 3H{pA +Pa) 
3H{pb +Pb) 



0, 
0. 



(12) 
(13) 



Equations ([7]), JH]), IH]), p2|) . and (fT3l) are not independent because of the Bianchi identities. The field equation ^ 
can be derived by using the other equations. 

III. CONDITIONS FOR THE AVOIDANCE OF GHOSTS AND LAPLACIAN INSTABILITIES IN THE 

PRESENCE OF TWO PERFECT FLUIDS 



In this section we study the stability of the flat FLRW background for the action ([6]) in the presence of two perfect 
fluids. 

Let us first consider scalar metric perturbations ^ , and $ with the line element (33 | 



(1 + 2'^)dt^ + 2^^x dtdx' + + 2$) dx^ 



(14) 



Here we have gauged away a scalar perturbation E that appears in the form E^ij in front of the term dx^. This fixes 
the spatial part of the gauge-transformation vector In the following we choose the uniform-field gauge for which 
the field cj) has no perturbations, that is, S(j> — 0. This fixes the temporal part of the vector <^^. 
The two perfect fiuids can be described in terms of the following Lagrangian 



5pf 



d'^Xy^ [pAiPAjSA) +Pb{PB,Sb)] , 



(15) 



where pi and Si {i = A, B) correspond to the chemical potential and the entropy per particle respectively. We will 
employ the method introduced in Ref. |35i] to study the perfect fluid from a Lagrangian point of view in order to 
extract the conditions for the absence of ghost and Laplace instabilities. In the following we will summarize and 
simplify the method given in Ref. [35j (see also Ref. (36j). 

Since we are interested in those fluids with equations of state of the kind pi = WiPi {i = A, B), we consider pi as 
functions of pi alone. Here the chemical potential of each fluid is defined as piUi = pi + Pi, where Ui — dpi/dpi is 
the number density of the fiuid i. In fact, it is sufficient to give the equations of state pi a n^'. We define the fluid 
4-velocity Ua associated with the chemical potential p, as Ua = p~^da(., where ^ is a scalar field. The normalization 



condition for the 4-velocity allows us to write the 4-velocity of the fluid A, as pA = y'—g'^PdaiAdpiA- After we 
perturb the field £a as iA{t) + SIa^ we can expand the matter Lagrangian at second order and then perform the field 
redefinition S£a — —PAVA, where va is chosen to represent the independent scalar degree of freedom of the fluid A. 
Along the same lines, the independent scalar degree of freedom for the fluid B corresponds to vb- Since, at linear 
order, the scalar fluids do not contribute to the vector perturbations, it is sufficient to study their scalar contributions 
alone in order to derive the conditions for the avoidance of ghosts and Laplacian instabilities. We will discuss this 
issue in Appendix A. 

We perturb the action S = J d'^x^—g C + S'pf up to the second order. After integrations by parts, the second-order 
action is given by 



5(2) = 



d^xa^ 



\ 2 



^ I + ^ (a$)2 - 3wi$2 

6 j 

pl{\ + Wl) 



-3$^p;(l -l-w;) [ill - iHwiVi) + -H^{l + wi)p, 



(16) 
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where I = A, B, and 

wi = 2iGi-2XGi^x)-2XiG5.x^H-G5,^), (17) 



W2 = -2Gz^xX(t) + AGiH -l&X^Gi.xxH + A{(t>Gi.^x -^HGi.x)X + 2G 



4, 



+ 8 X^HG^^^x + 2HX (6G5,^ - 5 Gs^x'i^H) - AG5,xx^X^H^ , (18) 

«;3 = iX{K^x + 2XK,xx) + &X{'iX4>HG:,,xx~G3,^xX~G^,^ + &H4>G:,,x) 

+ 18 (4 HX^G^^xxx - HG4 - 5 X0G4,0jf - G4,00 + 7 HG^.xX + 16 HX^G^.xx - 2 X^G4,4,xx) 
+6H^X{2 H^Gz^xxxX^ - 6 X^G^^^xx + UXH^G^^xx - 27G^^^xX + 15 iJ^Gg.x - 18G5^0) , (19) 

Wi = 2G4-2XG5,ci>-2XG5,x(t>- (20) 

The equation of motion for x gives rise to the following momentum constraint 

W2^' = 2iyi$+ pi{l + wi)vi, (21) 



l=A,B 



which can be used to integrate out the field ^E". After replacing in Eq. by using Eq. (PT|) . we need to integrate 
the terms i9^<i><i> and d^^vi by parts, so that the action reduces to 



where Qi — (<i>, va, vb), and Ay , Bij, Cij, Vij are the components of the 3x3 matrices A, B, C, V. 
Imposing the matrix A to be positive definite leads to the no-ghost conditions 

_ wi(4wiW3 + 9wi) 

= 5-^ >0, (23) 

and wa(1 + > 0, w_b(1 + wb)pb > 0. The latter two conditions are automatically satisfied for radiation 

{wa = 1/3) and non-relativistic matter (wb =0^). The speed of propagation Cs for the fields can be found by solving 
the following discriminant equation 

det {c^A - a^C) = . (24) 

This has two trivial solutions cj. = wa and cj. = wb, which are not negative for radiation and non-relativistic matter. 
In order to avoid the Laplacian instability associated with the remaining solution of Eq. (|24p we require that 

^2 ^ 3(2w^W2g - wlwj + 4wiW2Wi - 2w^W2) - 6w^[(l + wa)pa + (1 + wb)pb] ^ ^ ,25) 
~ wi(4-u;i-u;3 -f 9w|) ~ ' 

We also consider tensor perturbations /i^ characterized by Sgij — a^{t)hij, where hij is traceless (/i\ = 0) and 
divergence-free (/i*-' j — 0). We decompose hij into the two polarization modes, hij — ftgg e® -I- h^efj, where the 
matrices e®- and e® are normalized to be e^j{k)eij{—k)* = 2, (where A = ©,®) and efj{k)efj{—k)* = in Fourier 
space. The second-order action for tensor perturbations can be written as 



4'^ =E / d'^a^Qr 

\ ■J 



K-^ {dh 



a 



A J 



(26) 



A 

where the conditions for the avoidance of ghosts and Laplacian instabilities are given by 



Qt ^ ^>0, (27) 
4 ^ ^>0. (28) 

Wi 

For the consistency of the theories given by the action ^ the conditions (P^ . (|7f|) . and need to be 

satisfied. Finally, as we will show in Appendix A, the vector modes do not add any conditions to those derived above. 



5 



IV. APPLICATION TO EXTENDED GALILEON DARK ENERGY MODELS 



In this section we shall study the cosmology based on the extended Galileon models by taking into account the 
conditions for the avoidance of ghosts and Laplacian instabilities derived in Sec. IIIII 
The covariant Galileon without the field potential |TTI| corresponds to 

K=-c,X^ G3 = ^X, G4 = iMpV^^^ G5 = ^X^ (29) 

where Ci {i = 2,3,4,5) are dimensionless constants, Mpi is the reduced Planck mass, and M is a constant having 
the dimension of mass. In this case it is known that there exists a cosmological tracker solution characterized by the 
condition = constant fl^. Along the tracker the dark energy equation of state evolves as wue — —7/3 (radiation 
era) — > wde — ^2 (matter era) — > wde = — 1 (de Sitter era). However, the tracker solution for the covariant Galileon 
is not favored by the joint data analysis of Supernovae la (SNIa), Cosmic Microwave Background (CMB), and Baryon 
Acoustic Oscillations (BAO) [s^. This comes from the unusual evolution of wde during the matter era away from 
— 1. Only the solutions that approach the tracker at late times can be allowed observationally. 

Let us consider the generalization of the covariant Galileon to find tracker solutions with different wde. We take 
the following functions 

K = -C2M2^'''"^XP' , Gg = csMl-^P'XP' , G4 = -M^i - aM^-^P^XP^ , G5 = Sc^M-'^'+^P'^^ XP'^ , 

where and pi [i = 2,3,4,5) are dimensionless constants, and Mi [i = 2,3,4,5) are constants having dimensions 
of mass^. For two perfect fluids we consider radiation (energy density pA — Pr, equation of state wa — 1/3) and 
non-relativistic matter (energy density ps = Pm, equation of state wb — 0). For the choice (I30[) the field equations 
jl]) and ([8]) can be written as 

iH^Ml^ = pDE + + Pr , (31) 
(3if2 + 2H)Ml, = -pBE - Pr/3 , (32) 

where 

PDE ^ 2XK^x - K + eH^XGs^x - QH^Gi + 24H^X{Gi,x + XGi^xx) + 2H^X{5G5^x + 2XG5^xx) , (33) 
PUE ^ K - 24)XG3^x + 2(3ij2 + 2ij)G4 - UH^XGa^x - AHXGi^x - SHXGi^x - SHXXGi^xx 

-2X{2H^(j) + 2HH(j) + 3H'^4))G5,x - 4,H^4> X^Gs^xx , (34) 

and G4 = G4 - M2 /2 = -c^M^-^P'^XP'^. 

For the covariant Galileon {p2 = P3 — I, Pi — P5 — 2) there is a tracker solution with Hcj) — constant, in which case 
all the terms in Eq. (j33p are proportional to 0^. We search for tracker solutions characterized by the condition 

H,j>^'^ = constant , (35) 

where the power g is a real constant. If we choose the following powers, all the terms in Eq. p3p are proportional to 

P2=P, P3=P+{2q-l)/2, pi=p + 2q, p5=p+{6q-l)/2. (36) 

Note that the covariant Galileon corresponds to p = 1 and g = 1/2. 

Let us study whether the tracker really exists or not for the powers given by Eq. (f36| . Before doing so we first 
discuss conditions for the existence of a de Sitter solution characterized by </> = constant and H = constant. We 
introduce the following variable 

x^^. (37) 



^ Kimura and Yamamoto |26|| studied the model with p2 = 1,P3 ^ 0, C4 = 0, C5 = 0, which recovers the Dvali- Turner model |38|| . 
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For the simplification of the background equations it is convenient to define 

1/2 M. = MpL 



2p3 



where H^s is the Hubble parameter at the de Sitter solution. For the covariant Galileon one has M3 = M4 = A/5 — 
(i/|gA/pi)^/'^. From Eqs. (pijl and we find that the de Sitter fixed point where ^ = and = is present under 
the following conditions 



C2 - o — (3a " 4/3 + 2) , (39) 

^/2 / 2 X^^*^ 

C3 = 7 — [3(p + <?)(«- /3) + p] , (40) 



2p + g - 1 



"-dS 



where Xds is the value of x at the de Sitter solution, and 



<..i<^(*)"c„ ,.2Vi,J'W"\.. ,41, 



. 2 , 

In order to discuss the cosmological dynamics it is convenient to introduce the following variables 



a; \ 1 



p + 2q 
1 + 2, 



f]^ = — • (42) 



The solution psp . in terms of the variable ri, is given by ri = 1. The de Sitter fixed point corresponds to (ri, r2, il^) = 
(1,1,0). We will only consider cosmological dynamics for which ri and r2 are both positive at all times: otherwise the 
inverse relations of may be ill-defined. On the other hand, we will see later on that a viable cosmology cannot 
allow a change of sign for the variable r2. Defining the dark energy density parameter f^DE = PDE/(3iJ^Mp[), it 
follows that 



O ^FV2 

"DE = 



ri {n [I2(a ~ P){p + q) + ip - n (2p - l)(3a - 4/3 + 2)] - 3a(2p + 4g + 1)} + 4p{p + 3q + 1) .(43) 



The Friedmann equation gives the relation flm = Pm/i-^H^M^^) = 1 — fJ,. — fluE- For the initial conditions where ri 
is positive but small (0 < ri ^ 1), we require that 

^—^ > , (44) 

2<? + 1 - ' ' 

in order to have J^de — ^ 0. In the following we will replace rim with 1 — £7^ — J^de- 

We can obtain the differential equations for ri, r2, and f?^ after deriving (j) and H from Eqs. ([5]) and e.g., 
r[/ri = — 2g0/(iJ0) — H/H^, where a prime represents a derivative with respect to = Ina. The r.h.s. of these 
differential equations can be expressed in terms of ri, r2, and fir together with the coefficients a, /3, p, and q. 

The equation for ri is given by 

r'l = [nin - l){/3(2p(ri - l)(2ri + 1) - 6g(ri + 1)) + ri{3aip ~ pn + 2q) - 2pri)} 

x{2((3 + nr){2p - 1) + 12g)rf + 3(4/3(1 +p + 3q) - 3a(l + 2p + 4q)n 

+4((1 + 3a - 3/3)p + 3(a - ^)q)rl + (-2 - 3a + 4/3)(2p - I)r3)r2}]/A , (45) 

where 

A = 2[-2rf +^''"^^/'^+^*\3a(p + 2q)(2p + 4g - 1) - (2p + 2q - 1)((2 + 6a)p + 6aq)ri + (2 + 3a)p(2p - l)r?) 
_^4^^(i-p)/(i+2'?) (2p2 _ p _^ ^8^2 _^ ^2pq -3q- 3{2p + 2q-l){p + q)rj + 2p{2p - l)r?) 

+r?{9a2(2g(2gri -2q-l) + {2p^ + p + 8pq){ri - l))(ri - 1) + Ap^rf + 6apri{{2p + Aq + l)ri - 2p- 69- l)}r2 
+3/32(ri ~ l){4p(ri - 1)(1 + 2ri + 6q{ri + 1)^) + Ap'^in - 1)(1 + n{2 + Sn)) + 6g(ri + l)(2g(r^ - 3) - 2)}r2 
-3;3r^{4pri((ri - l)(2pri + 1) + 2q{rl - 2)) + a(ri - l){8p^{ri - l)(2ri + 1) + 6g(4gr? - (1 + 2q)ri - 3 - lOg) 
+ip{n - l)(12g(ri + 1) + n + 2))}r2] . (46) 
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In addition to Eq. (P5|) we write the remaining equations of motion for the variables r2 and $7^- We find that the 
second Einstein equation is equivalent to 

^{4(3g + + 2q) - 2r2{p + 2q)rf^ {3ari{p - q - 1) - 3^(p - 3g - 2) - rl[3ap - 3/3(p + q) + p + 3aq]}} 

+ '^{2q+ l){rf^r2{3/3ip + q) + rl[3ap - 3(3{p + q) + p + 3aq] - 3an{p + q)} + 2q} 

- 3{2q +l){p + 2q){rf^r2[rfi4l3 - 3a - 2) + 3an - 4^] + 2} - 2(2q + l){p + 2q) 17^ = , (47) 
whereas the equation of continuity for the radiation fluid leads to 

{2q +l)ip + 2q) ^ - 2(3<7 +l){p + 2q) ^ - 2q{2q + 1) ^ + 4(2g + + 2g) - . (48) 
s Ir ri r2 

These two equations, combined with Eq. (j45p . completely determine the dynamics of the system. 

A. Tracker 

From Eq. we find that there is a fixed point characterized by 

n = 1 , (49) 

which corresponds to the tracker solution with Hcj)^'' — constant. Considering the homogenous perturbations 5ri 
along the solution ri = 1, it follows that 

. , 6(p + 2g)-3 + (2p-l)a. + 3r2 ^ 
2(pr2 + 2q) 

For pr2 + 2g > the tracker is stable along the ri direction provided that 

6(p + 2g) + (2p-l)f7^>3(l-r2). (51) 

Along the tracker the dark energy density parameter is simply given by 

rioE = r2 . (52) 

For < ilDE < 1 the r.h.s. of Eq. ([?T|) is within the range < 3(1 — r2) < 3. As long as the condition ([?T|) is 
satisfied, which includes the case of the covariant Galileon, the solutions stay at the tracker. With the increase of r2 
the solutions finally reach the de Sitter fixed point characterized by ri = 1 and r2 = 1 . 
Along the tracker the equations for r2 and 17,. are given by 

jp + 2q){nr + 3 - 3r2) 

= ^^q ^''^ 

K - (54) 
pr2 + 2q 

Combining these equations, it follows that 

— -(l + s)7T^=4(l + s), where s=^. (55) 
r2 ilr 2q 

Integration of this equation leads to 

r2 = 010^(1+^)^,1+^ , (56) 

where ci is a constant. Since Qr oc a^^H^^, the evolution of the variable r2 is r2 oc i/^2(i+s)^ Since we want r2 to 
be subdominant at early times and to grow on the tracker solution, we further impose 

1 + 5 = 1 + ;^ >0. (57) 
2q 
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Substituting Eq. (|56p into we obtain the integrated solution 

ci a'*(^+') 17^^+" = 1-17^(1 -C2a), (58) 

where C2 is a constant. Therefore, the dynamics on the tracker depends only on the free parameter s = p/{2q) 
and the two initial conditions ci,2. For a particular choice of s such as s = 1 ( which corresponds to the covariant 
Galileon), it is possible to derive the explicit solution for fir in terms of a [itI. |37|. For general s, however, we cannot 
necessarily find an explicit expression for fir- The density parameter of dark energy along the tracker is given by 
^^DE = r2 = cia^'^^^'^^nr^^^ . Writing the density parameters of dark energy and radiation today (a = 1) as ilj^g and 
ni^-' respectively and using Eqs. ([5^ . ([5^ . and the coefficients ci and C2 are found to be 

Cl = uT^ , C2 = TTTT , (59) 

where fim'' = 1 — ~ ^r^^ density parameter of non-relativistic matter today. 

The dark energy equation of state wde = Pde/pde ^nd the effective total equation of state Wcs = ~1 ^ 2H / {SH^) 
along the tracker are 

3 + s{3 + nr) 3r2(.s + l)-fi, 

wde = ^7 — p— , Woff = — — — . 60) 

3(sr2 + 1) 3(sr2 + 1) 

In the early cosmological epoch in which the condition ftuE = r-2 <C 1 is satisfied, one has wde — — 1 — s(l + fir/ 3) and 
WeS — flr/3, respectively. Hence the evolution of wde during the radiation and matter eras is given by wde — — 1— 4s/3 
and Wde — ~1 — s, respectively. At the de Sitter fixed point (r2 — 1 and fir = 0) it follows that wde — WcS — ~1- 
The tracker solution for the covariant Galileon (s = 1) is incompatible with observations because wde is away from 
— 1 during the matter and radiation eras [37l] . For the compatibility with observations we require that 

s=f<l. (61) 
2q 

In the regime ^2 <C 1 the conditions for the avoidance of ghosts and Laplacian instabilities are given by 

^ ~ 6(?[p-3(«-2/3)g]r2>0, (62) 

4 ~ {^p^iflr + 3) - 2p'^{{flr + 3)(6/3 - 3a + 2) - 2q[3flr + H - 3(a - 2P){flr + 3)]} 

-3{l3{flr + 3) + Sq^iflr + 5)(a - 2/3) - 2q^{7flr + 27){a - 2/3) + q[3a{flr + 3) - 2/3(5^2^ + 17)]} 
-pii^r + 3){3a - 12/3 - 1) + Aq^[{a - 2f3){9flr + 33) - 2{flr + 5)] 

+g[12(2;3 - a){3flr + 10) + 6flr + 22]}}/[24g2(2p + Aq - l){p - 3{a - 2/3)q}] > , (63) 

§^ ^ i[2 + 3(a-2/3)r2]>0, (64) 

^2 ^ ^ 3{2[2(a - 2/3)g + 3/3]p + 8{a - 2/3)g^ + /3(16g - 3) + /3(2p + 4g - l)a.} ^ ^ 

"^■^ 4(j(2p + 4g-l) '"2 _ ■ l J 

These results reproduce those derived in Refs. (l6l. ITtI for the covariant Galileon {p — l,q — 1/2). Let us consider the 
case in which the parameters a, l3, p, q are not very different from the order of unity. Since r2 ^ 1 the conditions 
((64| and ((65| are automatically satisfied. We see here that r2 cannot change its sign, as this implies the violation of 
the condition Qs > 0. Since we only consider the case r2 > at the initial stage, the scalar ghost can be avoided for 

q[p- 3{a ~ 2/3)q] > . (66) 

For p and q satisfying q{2p + Aq — 1) > 0, the Laplacian instability of the scalar perturbation is absent as long as the 
numerator in Ecl (l63l) is positive. For the covariant Galileon this corresponds to the condition 8 + 10a — 9/3 + 51^(2 + 
3a - 3;3) > [ii, 13. 



B. de Sitter solutions (ri — r2 = 1) 



We study the stability of the late-time de Sitter solution characterized by ri = r2 = 1. At the de Sitter fixed 
point the system is effectively described by one scalar degree of freedom $, with the second-order action S'(^) = 
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/ (i^a;a^(5s[$^ — (c|/a^)(V$)^]. For homogenous perturbations (comoving wavenumber fc = 0) the scalar perturbation 
obeys the equation of motion -^{a^Qs^) — 0, whose solution is 



= ci + C2 / -jTT- dt , (67) 

where ci and C2 are constants. Now we are considering the de Sitter solution where (j) is constant, in which case Qs 
does not vary in time. Since a cx e^"*^*, the second term in Eq. (|67p decays in proportion to e~^^'^^^. Hence the de 
Sitter solution in our theory is classically stable against homogeneous perturbations. 

The conditions for the avoidance of ghosts and Laplacian instabilities (against inhomogeneous perturbations) at 
ri = r2 = 1 are given by 

05 ^ 6(p + 2g)(3a - 6/3 + 2) [p - 3(a - 2/3)g] 

M2, [2p - 6(a - 2/3)(j - 3a + 6/3 - 2]2 ^^^> 

4 = {6/3 + V + p [9(a - 2^)2 + 3a - 12/^ + 4g(6/3 - 3a + 2) - 2] + 3(a - 2/3)[3/3 + g(9a - 12/? - 8g + 6)]} 
3(2/3-a)(2g + l) + 2p-2 



> ('691 

6(6/3-3a-2)(p + 2q)(2p + 4g-l)(p-3a(7 + 6/3g) - ' ^ ^ 

^ = i(3a-6/3 + 2) >0, (70) 
2(2p + 4g- 1) -3a 



"pi 



(2p + 4g-l)(3a-6/? + 2) 



>0. (71) 



C. The solutions in the regime ri ^ 1 and r-2 <C 1 

Equation (|45ll shows that there is another fixed point characterized by 

ri=0. (72) 

Let us study the behavior of the solutions in the regime ri ^ 1 and r2 <C 1. In doing so we consider the parameter 
space with p > 1 and q > 0. Then the equations for ri, r2, and are approximately given by 

^ (3 + ».)(2p-l) + 12g 

- 2(2p + 6g-l) ^^^^ 
(p + 2g)[9-6p + r!,(7-2p+12g)] 
'^^ " 2(2g + l)(2p + 6,-l) ^'^^ 

n'^ ~ nrinr-i). (75) 

From Eq. (|75p there are two fixed points: fl^ = 1 and fJ^ = 0. During the radiation era {^Ir = 1) integration of 
Eqs. dZS]) and ^ gives 

4p + 6t;-2 (p + 2g)(S-4p + 6g) 

ri ~ a2p+6<j-i ^ ~ Q (2,+i)(2p+6q-i) ^ ^75j 

whereas during the matter era {Qr — 0) one has 

3(2p + 4i;-l) 3(p + 2q)(3-2p) 

ri ~ a2(2p + l39-l) ^ ~ £j2(2g+l)(2p + 6,-l) _ (77) 

The variable ri increases for p > 1 and q > 0, which is followed by the approach to the tracker solution. Whether 
the variable r2 grows or not depends on the values of p and q. If p > 3/2 and q > 0, for example, r2 decreases in the 
matter era. However the dynamics will in general change as the solutions approach the tracker, ri — >■ 1. 

The dark energy equation of state wbe and the total effective equation of state Wcs in the regime ri ^ 1 and 
r2 ^ 1 are approximately given by 

" 2(l-2p-6g) ' " S"" • 

If p = 1 and q = 1/2, we have wue = -(1 + ftr)/8 fl^flTt. 
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For {p — l)/{2q + 1) > the ghosts and Laplacian instabihties can be avoided for 

3{p + 3q){2p + 6q~ l)/?rJP-')/''«+'V2 > , 
p + 3q-2 



Ml, 



4 



Qt 



2{p + 3q){2p + 6q-l 



1 r 

4 



1 - 3;9r2ri' 



(p-l)/(29+l) 



>0, 



3(4p+12g-5-3»^) (p-i)/(2g+i) ^ 
4p + 129 - 2 ^ ~ ' 



(79) 
(80) 
(81) 
(82) 



In the regime ri <C 1 and r2 <C 1 the conditions ([5T|) and ([5^ are automatically satisfied. 

D. Other fixed points 

From Eq. (j45p we see that in general there are other two more complicated fixed points, ri = ra^b, those which 
satisfy the equation 



p{3a - 4/3 + 2)r2 + [2/3(p + 3q) - 3a{p + 2q)]rj + 2l3{p + 3q) = , 



(83) 



with j = a, b. Whether or not these fixed points are viable or not depends on their stability and the chosen parameters 
of the model. 

Then we will consider only those models for which either there are no real solutions to Eq. ([55]) . or, if they exist, 
they are placed outside the range of interest corresponding to the interval < ri < 1. Therefore we need to fulfill the 
following condition 



where 



ra,b 



(A < 0) V [A > A (r^ < V > 1) A (rb < V rfc > 1)] , 
3a{p + 2q) - 2(3{p + 3g) ± \/A 



2p (3a -4/3 + 2) 
A = [2/3(p + 3g)-3a(p + 2<7)]2-8/3p(3a-4/3 + 2)(p + 3g) 



(84) 

(85) 
(86) 



E. Theoretically viable parameter space 

Let us discuss the relevant constraints we have found so far in order to restrict the allowed parameter space. 
Therefore we have 



l + 2q 
P 



> 0. 



2q 
V_ 

2q 

p + 3q-2 



< 1: 

> 0. 



ip + 3q)i2p + 6q-l) 
{p + 3q){2p + 6q-l)P > 
3a - 6^ + 2 > 
{p + 2q)[p-3{a-2l3)q] >0 
q[p~3{a-2/3)q] >0 



41 



dc Sitter 



> 0. 



r \ > 

^Slri=l,r2<Cl,n^=0 - ■ 



order 


to 


have 


fiDESubdominant et 


early times, 


order 


to 


have 


a growing f^DE along 


the tracker 


order 


to 


have 


the tracker consistent with data, 


order 


to 


have 


'^•S' Li<Cl,r2<l - ^ ' 




order 


to 


have 


'3s' ri<l,r2<l ^ ^ ' 




order 


to 


have 


Qrldc Sitter > 0' 




order 


to 


have 


Qsldc Sitter > 0' 




order 


to 


have 


Qs\ri = l,r2<^l ^ ^ ' 





>o, 

r '^ \ > D 

''Slri = l,r2<Cl,Or = l — 



''Tide Sitter 
„2 



(87) 



(89) 

(90) 

(91) 
(92) 
(93) 
(94) 

(95) 
(96) 
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The conditions ([57])- ([101) impose 



g > ^ , and 1 < p < 2q . 

This imphes that negative values of p and q are excluded. Furthermore, the condition (j9ip implies 

/3 > 0. 

Altogether we find the following allowed parameter space 
1~ 



(97) 



(98) 



q > 



2 

2p-2 



A (1 < p < 2g) A 
P 



where 



bq + d 6q 



P 1. 



3a - 2p + 6aq + 2 



-2)A^</3</3„ 



6(2g+ 1) 



(99) 



(100) 



Although the above conditions have been derived for different initial conditions, we impose them to be true at the 
same time. This is because we do not know a priori initial conditions in the early Universe, as we have not specified 
how inflation works in these models. 

To be more concrete, let us consider the theory with p — 1 and q = 5/2. In this case the dark energy equation of 
state along the tracker is given by wde = — 1.2 during the matter era. The condition (|99|) reduces to 

(101) 



(a > 2/3) A (a > A ^ > A 15a < 2) V \^l5a >2Aa<2/3+ — A3a<22 

Adding the condition (|84l) . the parameter space reduces to 

(a > 2/3) A (a > A 2\/l7 ^3(272 - 75a)Q + 68 + 561/3 > 408a + 68 A 75a < 34 



V 75a > 34 A a < 2;S - 



15 



A 3a < 22 



(102) 



During the transition from the regime ri — 1 and r2 ^ 1 to the de Sitter attractor (ri = 1 and r2 = 1) the tensor 
propagation speed squared c|, can be negative, depending on the model parameters. For p = 1 and q ^ 5/2 we have 



110 + (22 - 15a - 99/3 - 33;3f7^)r2 + 3(33/3 - a)rl 
ll[2 + 3(a-2/3)r2](5 + r2) 



(103) 



along the tracker. Since the transition occurs at late times, the contribution of the radiation density parameter can 
be neglected in Eq. (|103l) . If |a| ^ 1 and |/3| ^ 1, then is close to 1. For positive a and /3 of the order of unity 
it happens that becomes negative. Taking the derivative of Eq. (jl03p with respect to r2, has a minimum (the 
second derivative (P{c'^)/dr2, at r2 = ''2,minj is always positive in the region < r2,niin < 1 a-nd when the conditions 
pOip are satisfied) at 



''2,1 



60a - 275/3 + 3VA~ 
55/3 - 12a - 594/32 ^ 297a/3 



where 



Arnin EE 13 [ISOa^ - 3a(197/3 + 8) + 22/3(21/3 + 5)] 



(104) 



(105) 



In Eq. (jl04p we have chosen the solution with r2 > 0, because another solution is always negative if the conditions 
(llOip are satisfied. 

Plugging the solution pTO)) into Eq. we find that 



36V55yS~ - 225a2 - 15a(69/3 + 20) + 66;3(45/3 - 23) + 44 
11(30;3- 15a + 2)2 



(106) 
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Figure 1: Allowed parameter space in the (a, /3) plane (the area in black) for p — 1 and q — 5/2. We have used all the available 
conditions considered so far, i.e. (|102p and p07p . 

Therefore the allowed parameter space is characterized by 

A„,i„ < V [A,ni„ > A (ra^min < V r2,min > 1)] V [A,„i„ > A (0 < r2,min < 1) A 41^^ > 0] . (107) 

The first inequality C01TL6S from setting T2 min 

to be an imaginary number such that there is no minimum for c^. The 
second inequality allows the minimum of to be outside the interested range for r2. 

In Fig. [T] we plot the allowed parameter space for p = I and q = 5/2 in the (a, (3) plane constrained by the conditions 

mnm and (fiuTH . 

F. Numerical simulations 

In order to confirm the analytic estimation given above, we solve the differential equations of ri, r2, and fir 
numerically for given values of a, /3, p, and q. In particular we study the cosmological evolution for p = I and 
q = 5/2, in which case the parameters a and (3 need to satisfy the constraints ()102|) and (|107|) . 

In Fig. [2] we plot the evolution of JIde, ^^m, ^r, Wcs, and wde versus the redshift z = 1/a — 1 for p = l,q = 5/2, 
a = 1, P — 0.45, and Xds = 1 with the initial conditions corresponding to the tracker from the beginning (ri = 1). 
For this choice of a and f3 the conditions (|102p and (|107p are satisfied. Note that the coefficients Ci {i — 2, - ■ ■ ,5) are 
known to be C2 = 9.60, C3 = 18.07, C4 = 4.36, and C5 = 7.20 from Eqs. ((39 |) - (|4T|) . Furthermore both the fixed points 
given in Eq. (|83p exist in the region ra^b > 1- Figure [2] shows that the sequence of radiation (woff = 1/3, fir = 1), 
matter (wes = 0, flm = 1), and de Sitter {wes = —1, fluE = 1) eras is in fact realized. The dark energy density 
parameter grows as floE = r2 (x = t^"^/^ toward the de Sitter attractor (ri = r2 = 1). As estimated by Eq. 

the evolution of wde during the radiation and matter eras are given by wde = —1.267 and wde — —1-2, respectively. 
After the end of the matter-dominated epoch the solution approaches the de Sitter attractor with wde = 

In Fig. [3] we show the variation of wue for the same model parameters as those used in Fig. [5] with a number 
of different initial conditions. The approach to the tracker occurs later for smaller initial values of ri. This can be 
clearly seen in the numerical simulations of the cases (a)-(c) in Fig.[3l The case (a) corresponds to the early tracking, 
whereas the case (c) to the late-time tracking with smaller initial values of ri . In the regime characterized by ri <C 1 
and r2 <C 1, the analytic estimation (|75|) gives wde = —0.0625 and Wde = —0.03125 during the radiation and the 
matter eras, respectively. These analytic values of wde are in good agreement with their numerical values for the 
late-time tracking solution (such as the case (c) in Fig. [3]). 

In Fig. [¥] the evolution of the quantity Qs is shown for the same model parameters and the initial conditions as 
those given in Fig. [31 For the tracker our numerical simulations show that Qs grows according to Eq. (j62p in the 
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Figure 2: Evolution of the density parameters JIde, fim, fir, Weff, and wdb versus the redshift z for p = 1, g = 5/2, a = 1, 
/3 = 0.45, and Xds = 1- The initial conditions are chosen to be r\ = 1, r2 = 10~^°, and ilr ~ 0.9998 a.t z = 1.76 x 10^. In this 
case the solution is on the tracker from the beginning. 




Figure 3: Variation of wde versus z for p = 1, q = 5/2, a = 1, (3 — 0.45, and a;ds = 1 with several different initial conditions. 
The solid line corresponds to the tracker with the initial conditions same as those given in Fig. (2] The initial conditions for the 
cases (a)-(c) are (a) n = 4.0 x 10~^ r2 = 5.0 x 10"^®, fir = 0.9998 at z = 1.82 x 10^ (b) n = 1.0 x 10"^ r2 = 1.0 x 10"", 
fir = 0.9998 at z = 1.76 x 10^, and (c) ri = 1.0 x 10"'^, r2 = 1.0 x 10"^, fir = 0.99995 at z = 6.64 x 10^, respectively. 



regime r2 ^ 1 (i.e. Qs cx r2 oc t^^^^), wliicli finally approaches the value ([55)) at the de Sitter solution. For the 
initial conditions with ri ^ 1 and r2 <C 1 we find that the early evolution of Qs is well described by Eq. ([79)) . i.e. 
Qs (X r2 (X ^19/32 during the radiation era. As we can see in Fig. HI the evolution of Qs shifts to that of the tracker 
after the solutions reach the regime around ri = 1. Provided that r2 > initially, Qs always remains to be positive. 
As we see in Eqs. ([M)) and ([5T|) the quantity Qt/M^^ is close to the value 1/4 in the regime r2 ^ 1, independently 
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Figure 4: Variation of Qs/Mpi versus z for the same model parameters and initial conditions as those given in Fig. [3] The solid 
line represents the tracker solution, whereas the cases (a), (b), and (c) correspond to the evolution for the initial conditions as 
those given in Fig. [31 
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Figure 5: Evolution of c| versus z for the same model parameters and initial conditions as those given in Fig. [S] The solid 
line represents the tracker solution, whereas the cases (a), (b), and (c) correspond to the evolution for the initial conditions as 
those given in Fig. |3l 

of the values of ri. Numerically we confirmed that, for both the initial conditions ri = 1 and ri ^ 1, Qt/M^^ stays 
the value around 1/4 until recently and then it finally approaches the value ([70)) at the de Sitter solution. Provided 
that 3q! — 6/? + 2 > the no-ghost condition for the tensor perturbation is always satisfied. 

In Fig. [5] we illustrate the evolution of c| for the same model parameters and the initial conditions as those given in 
Fig. El The analytic estimation for the tracker gives c|. = 0.639 and c|. = 0.515 during the radiation and matter 
dominated epochs respectively, which show good agreement with the numerical result in Fig.[S] Finally Cg approaches 
the value 9.48 x lO^'^ with a smooth transition from the matter era to the de Sitter epoch. In the regime ri <C 1 and 
r2 <C 1 the analytic estimation (|5n)) gives c| ~ (13/544) (1 + 17,.)- In fact the numerical simulations for the cases (a). 
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Figure 6: Evolution of Cy versus z along the tracker for (i) p = 1, g = 5/2, a = 1, /3 = 0.45, a;ds ~ 1, and (ii) p = 1, g = 5/2, 
a = 6.1, /? = 3, a;ds = 1- In both cases the initial conditions are chosen to be r\ = 1, r2 = 10""^'', and Q.r — 0.9998. In the case 
(ii) the condition (|107|) is violated, so that temporally becomes negative during the transition from the regime r2 1 to 
the regime r2 = 1. 

(b), and (c) reproduce this value before the solutions reach the tracker. Since < c|. < 1 from the radiation era to 
the de Sitter epoch, the Laplacian instability of the scalar perturbation is absent. 

Figure [S] depicts the evolution of (3^ for the tracker with two different combinations of a and /3. The case (i) 
corresponds to the model parameters and the initial conditions as those given in Fig. [5] As estimated by Eq. (|65p . the 
tensor propagation speed squared along the tracker is very close to 1 but slightly less than 1 in the early cosmological 
epoch. In the case (i) the model parameters satisfy the condition (|107p . so that takes a positive minimum value 
0.745 at r-i — 0.847. Then the solution finally approaches the de Sitter attractor with cip = 0.751. In the case (ii), in 
which the condition ()107p is violated, the minimum value of is negative. This shows that the condition ()107|) is in 
fact required to avoid the Laplacian instability of the tensor perturbation. We have also run our numerical code for 
the initial conditions with ri <C l,r2 ^ 1 and found that the evolution of is similar to that for the tracker. The 
only difference is that the tensor propagation speed is slightly superluminal in the regime ri <C 1,?'2 ^ 1- In fact, 
Eq. ((82|) shows that (3^ is slightly larger than 1 under the conditions (|97|) and (|98|) . 

We have also carried out the numerical simulations for other values of p, q and confirmed the accuracy of the 
analytic estimation. If we choose larger values of q for given p, the dark energy equation of state wde along the 
tracker approaches —1. In the limit where g — >■ cxd the scalar propagation speed squared (j69p at the de Sitter solution 
reduces to — ?> a — 2/3 and hence we require the condition a > 2/3. On the other hand the condition (j94| reduces to 
a < 2/3 for g cx). Hence, if wde along the tracker is close to —1. the allowed parameter space tends to be smaller. 

V. CONCLUSIONS 

In the Horndeski's most general scalar-tensor theories we derived conditions for the avoidance of ghosts and Lapla- 
cian instabilities associated with scalar, tensor, and vector perturbations. The four conditions p3|. ((25)) . (|27p . and 
([28| need to be satisfied for the theoretical consistency. Vector perturbations do not give rise to any additional 
conditions to those derived for scalar and tensor perturbations. 

The Horndeski's action covers most of the dark energy models proposed in literature and hence our formulas are 
general enough to apply them to concrete models with second-order field equations. We proposed new kinetically 
driven dark energy models described by the functions (|30)) , which cover the covariant Galileon as a specific case. For 
the choice of the powers -pi [i = 2, 3, 4, 5) given in Eq. ([36| we showed the existence of the tracker solution along which 
}i(\p-i = constant. Finally the solutions approach a stable de Sitter attractor at which = constant. 

Along the tracker the dark energy equation of state during the matter dominance is found to be iude = ^ 1 ^ p/ (2(7). 
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The covariant Galileon {p = 1 and q = 1/2) corresponds to wde = — 2 during the matter era, which is not favored 
from the combined data analysis of SNIa, CMB, and BAO. The extended Gahleon model we proposed can alleviate 
this problem because wde can be close to —1 for p smaller than q. 

We clarified the theoretically allowed parameter space in which the ghosts and Laplacian instabilities are absent. 
For p = 1 and q — 5/2 we carried out numerical simulations to check the evolution of the background quantities (like 
wuE and f^DE) as well as the quantities such as c|, c^, Qs, and Qt- As we estimated analytically, the dark energy 
equation of state for the tracker evolves as wde = —1-267 (radiation era), wde = — 1.2 (matter era), and wde = — 1 
(de Sitter era), see Fig. [3] For the initial conditions with ri <|; 1, wue starts to evolve from the value estimated by 
Eq. ([75]) . The approach to the tracker occurs later for smaller initial values of ri . 

For the values of a and /? which are inside the allowed parameter space, our numerical simulations show that c|, 
Cyi Qs, and Qt remain to be positive in the cosmic expansion history. Note that the condition (|107l) is important to 
avoid that becomes negative during the transition from the matter era to the de Sitter epoch. While we showed 
the cosmological evolution for one choice of p and q, we also confirmed that the analytic estimation is trustable for 
other values of p and q. In the limit that p/q ^ the dark energy equation of state for the tracker mimics that in 
the ACDM model. 

It will be of interest to see how the combined data analysis of SNIa, CMB, and BAO places constraints on the 
tracker solution in the extended Galileon models. In order to confront this model with the observations of large scale 
structure and weak lensing, we also need to study the evolution of matter density perturbations as well as gravitational 
potentials. We leave these issues for future work. 



ACKNOWLEDGEMENTS 



A. D. F. and S. T. were supported by the Grant-in- Aid for Scientific Research Fund of the JSPS Nos. 10271 and 
30318802. A. D. F. is also supported by the ThEP Center. S. T. also thanks financial support for the Grant-in- Aid 
for Scientific Research on Innovative Areas (No. 21111006). We thank Alexander Vikman for his initial collaboration 
on this work. S. T. is grateful to the organizers of sixth Aegean summer school and the Dark Universe Conference for 
their kind hospitalities during which a part of this work was done. 



Appendix A: Ghost conditions for the vector modes 



The r LO-g host conditions for the vector modes in the presence of two perfect fluids can be found by using the method 
in Ref. [35j- Let us consider the perturbed metric 

ds^ = -dt^ + a'^idt dx"- + a^{5ij + C,;j -I- Cj^i) dx'dx^ , (Al) 

where 7^.^ = = Ci^i. In order to describe the vector perturbations at linear order, for the perfect fluid we write 
fiu-^ — A'^OccBa, where A"^ and Ba are the velocity potentials for the fluid A. We can choose the background values 
for A'^ and Ba as follows: A^ = 0, and diBA = bf = arbitrary constant. Then, for a plane wave in Fourier space, we 
have Ciki = 0. By splitting B^ as B^ = 6f a;* + bfSBs, and choosing the gauge for which C± = = SB a, (where C± 
is the component of Ci perpendicular to bf), and the arbitrary background quantities such that bfki = 0, we find 
that bf is parallel to Ci and both are perpendicular to fc^. 

After expanding the action at second order in the fields and integrating out the auxiliary fields, we obtain 



S = Jd^x \QXAC^ + 2QX2bfC,5BB+Q22bfhf5B\^ , (A2) 
where Qn, Q12, and Q22 are time-dependent coefficients. The two no-ghost conditions can be written as 

y ^ ^5 W,{k/af[{l + WA)pA+PB{l + Ws)] ^ 

2 [2(1 + wa)pa + 2pb(1 + wb) + wi{k/af] 

4[2(1 + wa) PA + 2{1 + wb)pb + wi{k/aY] 

These conditions are satisfied for wi > (which corresponds to the condition (1271) ) and (1 -I- wa)pa > 0, and 
(1 +wb)pb > 0. In General Relativity with one single perfect fiuid, only the condition (jA3|) holds, and it agrees with 
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the result in Ref. [s^, when wi = Mpj. Hence the vector perturbations do not provide additional constraints to those 
derived for scalar and tensor perturbations. 
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